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Abstract. We investigate trions, paired states and quantum phase transitions in 
one-dimensional SU (3) attractive fermions in external fields by means of the Bethc 
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corrections for these physical quantities are presented in the strong attractive regime. 
Numerical solutions of the dressed energy equations allow us to examine how the 
different phase boundaries modify by varying the inter-component coupling throughout 
the whole attractive regime. The pure trionic phase existing in the strong coupling 
regime reduces smoothly by decreasing this coupling until the weak limit is reached. 
In this weak regime, a pure BCS-like paired phase can be sustained under certain 
nonlinear Zeeman splittings. 
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1. Introduction 

Recent experiments on ultracold atomic systems confined to one dimension (ID) P-H] 
have attracted renewed interest in Bethe ansatz integrable models of interacting bosons 
and multi-component fermions. The most recent experimental breakthrough is the 
realization of a ID spin-imbalanced Fermi gas of 6 Li atoms under the degenerate 
temperature [5]. This study demonstrates how ultracold atomic gases in ID may be used 
to create non-trivial new phases of matter, and also paves the way for direct observation 
and further study of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states [6j[7]. 

Three-component fermions exhibit a rich scenario, revealing more exotic phases 
[8HT5]. Notably, strongly attractive three-component ultracold atomic fermions can 
form three-body bound states called trions. Consequently, a phase transition between 
Bardeen-Cooper-Schrieffer (BCS) like pairing superfluid and trionic states is expected 
to occur in the strong attractive regime p ^ lT3 |[T6H22] . So far, most of the theoretical 
analysis has focused on the attractive strong coupling limit. A pertinent discussion in 
this context is what happens at intermediate and weak attractive coupling regimes 

In this paper we consider one-dimensional three-component ultracold fermions 
with 5-function interaction in external magnetic fields. From a mathematical point 
of view, this model was solved long ago by Sutherland [23] and Takahashi [21] through 
Bethe ansatz techniques. Recently, integrable models of three-component interacting 
fermions [T6 | [TT t l2"5] have received renewed interest in connection with ultracold atomic 
gases. Advanced experimental techniques newly developed allow to explore three- 
component Fermi gas with different phases of trions, dimers and free atoms [26H29] . 
Remarkably, direct observation of a trimer state consisting of fermionic 6 Li atoms in 
the three energetically lowest substates has been just reported in [30]. This opens up 
an opportunity to experimentally study such novel quantum phases of trions and pairs 
in ID three-component Fermi gases, providing a physical ground and stimulus for the 
investigation of their theoretical aspects. 

Our aim here is to expand on the theoretic knowledge of ID integrable model of 
three-component fermions by undertaking a detailed analysis of how the different phases 
modify as the inter-component interaction decreases, ranging from a strong to a weak 
regime. We obtain analytical expressions for the critical fields and construct the full 
phase diagrams in the weak coupling regime by solving the Bethe ansatz equations. We 
extend previous work on this model [16j|T9] to derive higher order corrections for these 
physical quantities in the strong coupling regime. Numerical solutions of the dressed 
energy equations show that the pure trionic phase existing in the strong coupling regime 
reduces as the coupling decreases and it disappears in the presence of external fields 
when the weak regime is approached. In contrast to the two-component interacting 
fermions [31], nonlinear Zeeman splitting may sustain a BCS-like paired phase in the 
three-component attractive fermions for weak coupling regime. 
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2. The model 

We consider a 5-function (contact potential) interacting system of N fermions with equal 
mass m, which can occupy three possible hyperfme levels (|1), |2) and |3)) with particle 
number N 1 , N 2 and A^ 3 , respectively. They are constrained to a line of length L with 
periodic boundary conditions. The Hamiltonian reads [23] 

h 2 N d 2 

H = - E -pOi + 9m E 6 ( x i ~ x i) + E z (1) 

i=l UX i l<i<j<N 

The first and second terms correspond to the kinetic energy and (^-interaction potential, 
respectively. The last term denotes the Zeeman energy E z = J2f=i ^e^ijig, 
with the Zeeman energy levels e l z determined by the magnetic moments fi B and the 
magnetic field B. For later convenience, the Zeeman energy term can also be written as 
Ez = —H 1 (N 1 —N 2 )—H2(N 2 —N 3 )+Ne, where the unequally spaced Zeeman splitting in 
three hyperfme levels can be specified by two independent parameters Hi = e— e l z (n B , B) 
and H 2 = e 3 z (fi B , B) — e, with e = J2a=i e z{^Bi the average Zeeman energy. 

The spin-independent contact interaction gin remains between fermions with 
different hyperfme states and preserves the spins in each hyperfme states, i.e., the 
number of fermions in each spin states are conserved. Although these conditions seem 
rather restrictive, it is possible to tune scattering lengths between atoms in different low 
sublevels to form nearly SU(3) degeneracy Fermi gases via broad Feshbach resonances. 
PQ[26j[27ll32]. Consequently the model still captures the essential physics relevant in 
the discussion of multiple phases in three-component ultracold Fermi gases. The inter- 
component coupling gin = —h 2 c/m with interaction strength c = — 2 /am is determined 
by the effective ID scattering length am [33J. It is attractive for gm < and repulsive 
for gm > 0. Here we will focus in the attractive case. For simplicity, we choose the 
dimensionless units of h = 2m = 1 and use the dimensionless coupling constant 7 = c/n 
with linear density n = N/L. 

The Hamiltonian (TjQ) exhibits spin £77(3) symmetry and was solved long ago by 
means of the nested Bethe ansatz [231 El]- In this approach its spin content was 
incorporated via the symmetry of the wavefunction. The energy eigenspectrum is given 
in terms of the quasimomenta {fc;} of the fermions through 

N 

£ = E *?. (2) 

3=1 

satisfying the following set of coupled Bethe ansatz equations (BAE) [2"3"]|2"3] 

it A a -k e + ic/2 = ^ A a -A^ + ic || A a - X e -ic/2 
JLi K-ke-i c/2 ^ A a — A^ — i c/ = \ A a - X e + i c/2 ' 

M a m - A £ - i c/2 l\ A M - Xt - ic' U 
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written also in terms of the rapidities for the internal hyperfine spin degrees of freedom 
A Q and A M . Above j = 1, . . . , N, a = 1, . . . , M 1; fi = 1, . . . , M 2 , with quantum numbers 



3^3 2^2]^! 



a three- 



Mi = A" 2 + 2 A3 and M 2 = N 3 . For the irreducible representation 
column Young tableau encodes the numbers of unpaired fermions (Ni = A^ 1 — A^ 2 ), 
bound pairs (A" 2 = N 2 — N 3 ) and trions (N 3 = A^ 3 ). 

3. Ground states 

For attractive interaction, the BAE allow charge bound states and spin strings. In 
particular, the SU(3) symmetry carries two kinds of charge bound states: trions and 
pairs. In principle, different numbers of unpaired fermions, pairs and trions can be 
chosen to populate the ground state by carefully tuning Hi and H 2 . 

In the strong coupling regime L\c\ 3> 1, the imaginary parts of the bound states 
become equal-spaced, i.e., a trionic state has the form {kj = Aj ± i|c|, A^} and for the 
bound pair {k r = A r ± i|c|/2}. Substituting these root patterns into the BAE Q, we 
find their real parts, from which the ground state energy in the strongly attractive 
regime can be obtained [16] 



E _ 7r 2 n? ( t 8n 2 + 4n 3 t 12 M i ^ 2 \ n 2 c 2 




L 3 V \ c \ 

ix 2 nl ( 12m + Qn 2 + 16n 3 1 , o 2 

+ —r- 1 H ^ V 7T^( 6n i + 3n 2 + 8n 3 ) - 2n 3 c 

6 \ 3|c| 3c z J 

ix 2 n\ ( 12m + 32n 2 + 18n 3 1 . n9 \ . , 

+ V i 1 + —^SR + 2W (fal + 16 " 2 + 9 " s) ) ■ (4) 

Here n a = N a /L (a = 1,2,3) is the density for unpaired fermions, pairs and trions, 
respectively. This state can be considered as a mixture of unpaired fermions, pairs and 
trionic fermions, behaving basically like particles with different statistical signatures [34J. 
For strong attractive interaction, trions are stable compared to the BCS-like pairing 
and unpaired states. From (T4|) we can obtain the binding energy for a trion, given by 
e t = h 2 c 2 /m and the pair binding energy, which is = h 2 c 2 /Am. 

In the weak coupling regime L\c\ 1 the imaginary parts iy of the charge bound 
states are the roots of Hermite polynomials Hk of degree k. Specifically, Hk(J^y) = 0, 
with k = 2, 3 for a bound pair and a trion, respectively [35]. The real parts of the 
quasimomenta deviate smoothly from the values evaluated for the c = case. With 
this root configuration, the ground state energy in the weak attractive regime can be 
obtained 

E re 2 , , , o N 

_~_( n 3 + 2n 3 + 3n 3 ) 

+ n 2 {ni{n 2 + n 3 )(ni + n 2 + n 3 ) + 2n 2 n 3 (n 2 + n 3 )) 

- 2\c\(nin 2 + 2ri\n 3 + 4n 2 n 3 + n 2 + 3nl). (5) 

The ground state energy (jS]) is dominated by the kinetic energy of composite 
particles and unpaired fermions and has an interaction energy consisting of density 
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density interaction between charge bound states and between charge bound states and 
unpaired fermions. From equation (JSJ) we can obtain the binding energy for a trion, 
given by St = 3h 2 \c\/mL and the pair binding energy, which is Sb = h 2 \c\/mL. For weak 
attractive interaction, the trionic state is unstable against thermal and spin fluctuations. 
This becomes apparent in the weak coupling phase diagrams presented in Fig. [TJ[d), Fig. 
E^c) and Fig. |3tc) below. 

4. Dressed energy formalism 

In the thermodynamic limit, i.e. L, N — > oo with N/L finite, the grand partition 
function Z = tr(e^ n ^ r ) = e~ G / T is given in terms of the Gibbs free energy G = 
E + Ez — (iN — TS, written in terms of the Zeeman energy E%, chemical potential 
/i, temperature T and entropy S [T6 1I36H38] . The Gibbs free energy can be expressed 
in terms of the densities of particles and holes for unpaired fermions, bound pairs and 
trions, as well as spin degrees of freedom, which are determined from the BAE ([3]). Thus 
the equilibrium state is established by minimizing the Gibbs free energy with respect to 
these densities. This procedure leads to a set of coupled nonlinear integral equations, 
from which the dressed energy equations are obtained in the limit T — > (see [T61I38] 
for details) 

e ( 3 )(A) = 3A 2 - 2c 2 - 3(i - a 2 * e«(A) 

-[a 1 + a 3 ]*e (2) (A)-[a 2 + a 4 ]*e (3) (A) 

e ( 2 )(A) = 2A 2 - 2(i - °- - H 2 - ai * e«(A) 

-a 2 *e 2 (A)-[a 1 + a 3 ]*e( 3 )(A) (6) 
= k 2 - (i - H x - a x * e^ 2 \k) - a 2 * e^A;). 

Here e- a \a = 1,2,3 are the dressed energies for unpaired fermions, bound pairs and 
trions, respectively and a,j(x) = 'hhjjcp^+x 1 ' s y m b°l "*" denotes the convolution 
dj * e( a \x) = f-Qa a j( x ~ y) e ^ a \y)dy with the integration boundaries Q a given by 
e^(±Q a ) = 0. The Gibbs free energy per unit length at zero temperature can be 
written in terms of the dressed energies as G = J2a=i f~ $-Qa ^(x)dx. 

The dressed energy equations (jSJ) can be analytically solved just in some special 
limits. In particular, they were solved in [TH] for strongly attractive interaction through 
a lengthy iteration method. Here we numerically solve these equations to determine the 
full phase diagram of the model for any value of the coupling. This allows to examine 
how the different phase boundaries deform by varying the coupling from strong to weak 
regime. The numerical solution is also employed to confirm the analytical expressions 
for the physical quantities and the resulting phase diagrams of the model in the weak 
coupling limit (see [19] for a similar discussion in the strong regime). 
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5. Full phase diagrams 

Basically, there are two possible Bethe ansatz schemes to construct the phase diagram 
of the system. One possibility is to handle with the dressed energy equations (jSJ). This 
approach was discussed for the strong attractive regime in [16J, where expressions for 
the fields in terms of the densities were obtained up to order of l/|c|. Alternatively, one 
can handle directly with its discrete version (equations (|2]) and (J3])) by solving the BAE. 
We adopt this second strategy here. In order to obtain the explicit forms for the fields 
in terms of the polarizations we consider the energy for arbitrary population imbalances 



E/L = iin + G/L + n l H l + n 2 H 2 , (7) 

which coincides with the ground state energy (T5]) obtained by solving the BAE ([3]). Then 
the fields Hi and H 2 are determined through the relations 

H x - 9E/L , E t - 9E/L (8) 
dni ' dn 2 

together with the constraint 

n = rii + 2n 2 + 3n 3 . (9) 
In the strong coupling regime, using the ground state energy (TJ]) we find 

00 f Arii 8n 2 An? 12, n0 8ni , \ 

1 = 1 ( 1 "9|^ + ^f + ^f + ^ (2n2 + ^ ~~ 3? (2 ™ 2 + ^ ) 

n 2 n 2 3 ( ^ Ani ^ 32n 2 ^ An 3 ^ (Qrii + lQn 2 + 9n 3 ) 2 2n 3 (Qrii + lQn 2 + 9n 3 )^ 



3|c| 9|c| 3|c| 27c 2 9c 2 J 
10n 2 nl ( 6ni + 3n 2 + 8n 3 \ 2c 2 

+ — -) + -- 

n 2 n 2 ( An x 40n 2 16n 3 (6ni + 3n 2 + 8n 3 ) 2 lAn^m + 3n 2 + 8n 3 ) \ 

H2 - — { l + ^\ + W\ + W + 3^ 27^ ) 

2ix 2 n\ An x + 32n 2 + 8n 3 t (6m + 16n 2 + 9n 3 ) 2 10n 3 (Qm + 16n 2 + 9n 3 ) 



9 V 3|c| 9|c| 9\c\ 27c 2 27c 2 J 

+ 16n 2 nl ( + 6(2n 2 + n 3 ) \ + 5c 2 
9|c| \ |c| / 6 



These equations provide higher order corrections to those derived in [16] using the 
dressed energy equations. For determining the full phase boundaries, we also need the 
energy-field transfer relation between the paired and unpaired phases Hi — H 2 /2, which 
can be extracted from the underlying two-component system with SU(2) symmetry 
[3T|l3lj. These equations determine the full phase diagram and the critical fields 
activated by the fields Hi and H 2 . 

Fig. Ufa) shows the ground state energy versus Zeeman splitting parameters Hi 
and H 2 determined from equation (jlj) with the densities m and n 2 obtained from (fTUj) . 
There are three pure phases: an unpaired phase A, a pairing phase B and a trion 
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Figure 1. Ground state energy vs Zccman splitting for different coupling values (a) 
strong interaction I7I = 10, (b) |7| = 5, (c) I7I = 1 and (d) weak interaction |-y | = 0.5. 
The white dots correspond to the numerical solutions of the dressed energy equations 
©. The black lines in (a) and (d) are plotted from the analytical results (flDj) and 
(fTTj). respectively. A good agreement is found between the analytical results and the 
numerical solutions in the strong and weak regimes. The pure trionic phase C, present 
in the strong coupling regime, reduces smoothly as I7I decreases and is suppressed in 
the weak limit. 



phase C and four different mixtures of these states. For small H±, a transition from a 
trionic state into a mixture of trions and pairs occurs as H 2 exceeds the lower critical 
value iff 1 . When H 2 is greater than the upper critical value if| 2 , a pure pairing phase 
takes place. Trions and BCS-like pairs coexist when H 2 l < H 2 < H 2 2 . These critical 
fields, derived from equation (ITU1) . are given by H^ 1 m n 2 llf(l + 2WI ~ W^)) 
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and Hf n 2 (^f- + ^(1 + ^ - 3^2 )) • The phase transitions from B ^ A + B ^ A 
induced by increasing Hi are reminiscent of those in the two-component systems [3l"tl4"T] . 
Basically in this region the highest level is far away from the other two levels, so 
the system reduces to the spin- 1/2 fermion case. The mixed phase containing BCS- 
like pairs and unpaired fermions can be called a FFLO phase. We mention that a 
discussion about the pairing nature of ID many-body systems can be found, for instance, 
in [391110]. For small H 2 , a phase transition from a trionic into a mixture of trions 
and unpaired fermions occur. Using equation ( flOj) . we find that the trionic state with 
zero polarization n\/n = forms the ground state when the field H < Hf~, where 
Hf 1 m n 2 (^j- — |^(1 + + g^r))- When Hi is greater than the upper critical value 

H± 2 n 2 + 7r 2 (l — gpbOJj all trions are broken and the state becomes a normal 
Fermi liquid. 

At intermediate coupling regimes, it is not possible to construct the full phase 
diagrams analytically. However, they can be determined by numerically solving the 
dressed energy equations (jSJ), as illustrated in Figs. Hb),[Tfc); |2fa), [2](b) andEJa), Efb), 
for the intermediary values of the coupling |c| = 5 and |c| = 1, respectively. The different 
phase boundaries modify slightly by varying the inter-component coupling through the 
whole attractive regime. In particular, the pure trionic phase existing in the strong 
coupling regime reduces smoothly by decreasing this coupling until it is completely 
suppressed. A careful numerical analysis of the phase diagrams for n = 1 and different 
values of |c| between |c| = 1 and |c| = 0.5 indicates that the critical coupling value at 
which the trionic phase disappears is around c c ~ 0.6. Other mixed phases involving 
trions, specially the phase (B + C) also reduce by decreasing |c|. 

In the weak coupling regime we obtain the expressions between the fields and the 
polarizations using equations (JSJ), (jHl) and Q 



Hi = —{2n\ + n\ + 4n x n 2 + Anin 3 + 2n 2 n 3 ) 
o 

+ ^(2n! + n 2 ). 

7T 2 

H 2 = — {n\ + 2n\ + 2n x n 2 + 2riin 3 + 4n 2 n 3 ) 

+ ^(2n 2 + m). (11) 

These equations together with the energy- field transfer relation Hi — H 2 /2 determine 
the full phase diagram and the critical fields activated by the Zeeman splitting Hi and 
H 2 . We observe that the density of trions n 3 does no appear independently in eqs. lfTTTl . 
in contrast to the corresponding equations in the strong regime fflQj) . Fig. QJd) presents 
the ground state energy versus the fields Hi, H 2 while Figs. [2](c) and [3|c) show the 
polarizations nx/n and n 2 /n in terms of Zeeman splitting, respectively. Now in the 
weak coupling regime there are just six different phases in the Hi-H 2 plane: we observe 
the disappearance of the pure trionic phase C in the presence of the fields, i.e., the trionic 
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(a) (b) (c) 



Figure 2. Phase diagram showing the polarization ni/n versus the fields Hi and H2 
for different coupling values (a) I7I = 5, (b) |7| = 1 and (c) weak interaction I7I = 0.5. 
The white dots correspond to the numerical solutions of the dressed energy equations 
(|6]). A good agreement is found between the analytical results (fTTj) (black lines) and 
the numerical solution in the weak regime. At intermediary coupling regimes, the pure 
trionic phase reduces by decreasing the coupling and it is not present in the weak 
regime. A zoom around the origin is presented in Fig. 2(b) (Fig. 2 (c)) to show the 
presence (absence) of the trionic phase. 

state is unstable against thermal and spin fluctuations. This behaviour is in contrast to 
the strong coupling regime, where the phase C is robust and trion states populate the 
ground state for a considerable interval of the fields. In addition, the phase where trions 
and pairs coexist (B + C) reduces significantly compared to the strong coupling regime. 
Interestingly, in contrast to the weak attractive spin-1/2 fermion system, a pure paired 
phase can be sustained under certain Zeeman splittings. For certain tuning Hi and H2, 
the two lowest levels are almost degenerate. Therefore, the paired phase naturally occurs 
and is stable. The persistence of this phase is relevant for the investigation of phase 
transition between BCS-like pairs and FFLO states. All these boundary modifications 
occur smoothly, as shown by a numerical analysis of the phase diagrams for different 
values of the coupling across all regimes. This indicates that all phase transitions in the 
vicinities of critical points are second order. This conclusion is consistent with previous 
analytical results [311 E]. We also mention that quantum phase transitions between 
different superfluid phases have been discussed in [12"] . 

We perform a similar analysis as in the previous strong case to extract the critical 
fields. Since the trionic phase C disappears for non- vanishing fields, less critical fields are 
found compared to the strong coupling case. For small Hi a transition from a mixture 
of trions and pairs into a pure paired phase occurs as H 2 exceeds the critical value 
H% ~ n 2 + . The transition from a mixture of trions and unpaired fermions into 

a normal Fermi liquid phase occurs as Hi exceeds the critical value H^ fa n 2 C^- + . 
The phase transitions B — > A + B — > A are reminiscent of those in spin-1/2 fermion 
systems. However, in the weak attractive two-component case the pure BCS-like paired 
phase is suppressed [31] and consequently it is not possible to investigate the phase 
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Figure 3. Phase diagram showing the polarization riijn versus the fields H\ and H2 
for different coupling values (a) I7I = 5, (b) |7| = 1 and (c) weak interaction I7I = 0.5. 
The white dots correspond to the numerical solutions of the dressed energy equations 
([6]). The black lines plotted from the analytical results (jlll) are in good agreement 
with the numerical results in the weak regime. At intermediary coupling regimes, the 
pure trionic phase reduces by decreasing the coupling and it is not present in the weak 
regime. A zoom around the origin is presented in Fig. 3(b) (Fig. 3(c)) to show the 
presence (absence) of the trionic phase. 

separation between a BCS-like paired phase and a FFLO state, in contrast to the three- 
component case, where this study is still possible in a weak regime. 

6. Conclusion 

We have studied the three-component attractive ID Fermi gas in external fields through 
the Bethe ansatz formalism. New results for the critical fields and complete zero 
temperature phase diagrams have been presented for the weak coupling regime. Previous 
work on this model has been extended to derive higher order corrections to these physical 
quantities in the strong regime. We have further confirmed that the system exhibits 
exotic phases of trions, bound pairs, a normal Fermi liquid and mixture of these phases 
in the strongly attractive limit. We have also shown how the different phase boundaries 
deform by varying the inter-component coupling across the whole attractive regime. In 
particular, the trionic phase that may occur in the strong coupling regime for certain 
values of the Zeeman splittings reduces smoothly by decreasing the coupling until the 
weak limit is approached, when the trionic phase is suppressed. Interestingly, in the 
weak regime, a pure paired phase can be maintained under certain nonlinear Zeeman 
splittings, in contrast to the two-component attractive ID Fermi gas. Our high precision 
of critical phase boundaries pave the way to further investigate quantum criticality in 
three-component interacting Fermi gas through the finite temperature Bethe ansatz. 
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